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Abstract. We present and discuss the high energy (E >4 keV) XMM-Newton spectrum of the Seyfert 2 galaxy, 
NGC 1068. Possible evidence for flux variability in both the neutral and ionized reflectors with respect to a 
BeppoSAX observation taken 3.5 years before is found. Several Fe and Ni emission lines, from both neutral and 
highly ionized material, are detected. The intensity of the iron Kct Compton shoulder implies that the neutral 
reflector is Compton-thick, likely the visible inner wall of the Nh > 10 25 cm -2 absorber. From the equivalent 
width of the ionized iron lines a column density of a fewxlO 21 cm -2 is deduced for the hot ionized reflector. 
Finally, an iron (nickel) overabundance, when compared to solar values, of about 2 (4) with respect to lower Z 
elements, is found. 
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1. Introduction 

The archetypal Seyfert 2 galaxy, NGC 1068, has been ex- 
tensively studied at all wavelengths. It was the discovery 
of broad lines in the polarized flux of this source that 
led Antonucci & Miller (1985) to propose the Unification 
model for Seyfert galaxies, which is at present the most 
popular scenario for this class of sources. In X-rays, af- 
ter the pioneering Einstein and EXOSAT observations 
(Monier & Halpern 1987; Elvis & Lawrence 1988), the 
GINGA discovery of a strong iron line (Koyama et al. 
1989), being interpreted as due to the reprocessing by cir- 
cumnuclear matter of the otherwise invisible nuclear radi- 
ation, brilliantly confirmed the Antonucci & Miller model. 
ASCA (Ueno et al. 1994; Iwasawa et al. 1997; Bianchi et 
al. 2001) demonstrated that the line is complex, implying 
reflection from both neutral and ionized material. The two 
reflectors model was confirmed, observing the hard X-ray 
continuum, by BeppoSAX (Matt et al. 1997a), which was 
also able to put a lower limit of about 10 25 cm -2 to the 
column density of the absorbing matter. Guainazzi et al. 
(1999) and Bianchi et al. (2001) demonstrated that the 
situation is even more complex, the line spectrum requir- 
ing at least three reflectors, one neutral (cold), one mildly 
ionized (warm) and one highly ionized (hot). Comparing 
the two BeppoSAX observations, taken about one year 
apart, Guainazzi et al. (2000) revealed a flux variability, 
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which they interpreted as due to a variation of the inten- 
sity of the hot ionized reflector, so placing an upper limit 
to its size of the order of a parsec or so. 

High spatial and spectral resolution Chandra and 
XMM-iVeiufon observations have confirmed the complex- 
ity of the circumnuclear region. The Chandra image 
(Young et al. 2001) revealed, especially in soft X-rays, 
a very rich morphology, with the brightest spot, how- 
ever, confined in a 1".5 (corresponding to 118 pc assum- 
ing H o =70 km/s/Mpc, Young et al. 2001) region around 
the nucleus. Indeed, the iron emission as well as the neu- 
tral and ionized continuum reflection are mostly confined 
in the nuclear region (Ogle et al. 2003). Many off-center 
point-like sources are also detected (Smith & Wilson 
2003). 

Gratings observations (Kinkhabwala et al. 2002; Ogle 
et al. 2003) have measured a line spectrum which is 
consistent with photoionized plasma with a wide range 
of ionization parameters. The line fluxes require signifi- 
cant contribution from resonant scattering (Kinkhabwala 
et al. 2002), as predicted by Band et al. (1990) and 
Matt et al. (1996). Given the limited energy range of the 
XMM-iVeuifon/RGS and the modest effective area of the 
Chandra /HETG at high energies, these results mainly 
concern the soft X-ray line spectrum. 

In this paper we analyse and discuss the high energy 
(>4 keV) XMM-Newton spectrum of NGC 1068. While 
the CCDs of the EPIC cameras have a significantly poorer 
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Fig. 1. Left: The whole band Chandra ACIS image of the nuclear region of NGC 1068. Right: The same, but only 
above 4 keV. Note that the extended emission is less prominent, and that one off-center source, CXJ024238. 9-000055, 
becomes by far the most prominent one. 



energy resolution than the Chandra HETG (the only grat- 
ing instrument currently working at the iron Ka energy), 
for many purposes the much larger effective area above 4 
keV overcompensate for the lower resolution. As an ex- 
ample, one can compare the wealth of informations on 
Fe and Ni K lines derived from the XMM-Newton ob- 
servation of the Circinus Galaxy (Molcndi et al. 2003) 
with the much poorer informations on the same lines (in- 
deed many of them not even detected) derived from the 
Chandra /HETG spectrum (Sambruna et al. 2001). 

2. Data reduction 

2.1. XMM-Newton 

NGC 1068 was observed twice by XMM-Newton on July 
29, 2000 and July 30, 2000 with the imaging CCD cam- 
eras, the EPIC-pn and MOS, adopting the Medium fil- 
ter and operating in Large Window (pn), Full Window 
(MOS1) and Prime Partial W2 (MOS2). We retrieved the 
ODF files from the public archive, and reduced them with 
the latest available SAS version (5.4.1), therefore using 
the most updated calibrations. The observed count rate is 
lower than the maximum defined for a 1% pileup for the 
selected pn subarray (see Table 3 of XMM-Newton Users' 
Handbook), but it is close to that value, so we conserva- 
tively decided to exclude patterns higher than in data 
reduction. An analysis with the epatplot tool indicates 
that same residual pile-up effect may still be present, so 
we extracted a spectrum removing the inner 6" circle re- 
gion. No significant differences are found, however, in the 
continuum and in the lines alike, if the comparison is lim- 



ited to the energy range of interest here, i.e. above 4 keV. 
We therefore concluded that, for the pn, pile-up is not 
relevant for our purposes. On the other hand, we shall 
not use MOS data because MOS1 has been used in Full 
Window mode and therefore strongly suffers from pile-up 
problems, while MOS2 does not add much information, 
expecially because we will limit our analysis to E >4 keV. 
Moreover, the pn appears to be the best calibrated detec- 
tor at the energies under consideration, as deduced from 
the analysis of the iron line complex in the spectrum of 
the Circinus Galaxy (Molcndi et al. 2003). 

Since the two observations are nearly contiguous and 
no significant flux variations are found in their lightcurve, 
we verified that the two spectra were consistent with each 
other and decided to combine them. We used the SAS tool 
merge to obtain a single event file and then extracted 
a spectrum from a region of 40" radii. The total expo- 
sure time, after rejecting time intervals of flaring particle 
background, is 61 ks. Spectra were analysed with Xspec 
11.2.0. 

2.2. Chandra 

Given the complexity of the morphology of the circumnu- 
clear region of NGC 1068, we also analysed the Chandra 
observation of NGC 1068 performed on February 21, 2000 
with the Advanced CCD Imaging Spectrometer (ACIS- 
S), for a total exposure time of 48 ks. The default 3.2 
s frame time, though resulting in the nuclear spectrum 
spoiled by pileup, does not affect our qualitative imag- 
ing analysis and the spectrum of the much fainter source 
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Fig. 2. Left: The whole band XMM-Newton EPIC-pn image of the nuclear region of NGC 1068. Right: The same, 
but only above 4 keV 



CXJ024238.9-000055 (see below). Data were reduced with 
CIAO 2.3, using CALDB 2.21. Spectra were analysed with 
Xspec 11.2.0. 

3. Spatial analysis 

In Fig.nthe Chandra image over the whole working band 
(left panel) and above 4 keV only (right panel) is shown. 
The detailed spatial analysis is discussed by Young et al. 
(2001), while the X-ray source population by Smith & 
Wilson (2003). The XMM-Newton image (whole band, 
left panel, and above 4 keV, right panel) is shown in 
Fig-Hi I n a U images, the r=40" circle over wich the XMM- 
Newton spectrum has been extracted is shown. 

For our own purposes, it is important to note that: a) 
the nucleus is partly extended over several arcseconds even 
above 4 keV, the extended emission accounting for about 
a quarter of counts in this band. Our 40" XMM-Newton 
extraction regions encompasses this extended emission, b) 
At a distance of about 28", a bright point-like source 
(CXJ024238.9-000055) is clearly apparent in the Chandra 
image (hereinafter we refer to it for simplicity as X— 1); in 
the hard band, it is by far the brightest off-nuclear point- 
like source. Smith & Wilson (2003) derived for this source 
a 0.4-8 keV flux of about 5xl0 -13 erg cnr 2 s^ 1 , corre- 
sponding to a luminosity, at the distance of NGC 1068, 
of 1.4xl0 40 erg s _1 . The spectrum (which we reanal- 
ysed) is quite hard, being consistent with either a very 
flat (r=0.85) power law or with multicolor disc blackbody 
with an inner temperature of about 4 keV. While we de- 
fer the reader to Smith & Wilson (2003) for a discussion 
of its possible nature, we note here that this source, be- 



ing so hard, may significantly contaminate the spectrum 
of the NGC 1068 nucleus in the energy range under con- 
sideration in this paper. Indeed, while the source is not 
readily visible in the whole band XMM-Newton image, 
it is clearly present in the E >4 keV image. To evaluate 
its contribution to the overall spectrum, we extracted a 
20" spectrum centered on NGC 1068 and a 10" spectrum 
centered on X-l. The flux we derived for NGC 1068 is 
14 times larger than that of X-l (which in turn is about 
half that measured by Chandra, Smith & Wilson 2003), 
and this is certainly a lower limit as the spectrum of X- 
1 is heavily contaminated by that of NGC 1068. In fact, 
in the X-l spectrum a broad iron line is clearly present, 
well fitted by two narrow components, one at 6.4 keV 
(neutral iron) and one at 6.7 keV (He-like iron). If we 
assume that they are the result of the contamination by 
NGC 1068 (no iron lines are apparent in the, admittedly 
poor, Chandra spectrum), this can provide a method to 
estimate the level of contamination of X-l by the nucleus. 
The EWs for these two lines are about 2.5 smaller than in 
the NGC 1068 spectrum (when calculated with respect to 
the whole continuum; see next section). Therefore, about 
40% of the spectrum extracted around X-l is actually due 
to NGC 1068, which implies that X-l does not contami- 
nate the 40" spectrum of NGC 1068 by more than about 
5%. Considering also that we are mostly interested in the 
emission lines, which are unlikely to be produced by X-l, 
we decided that this contamination is not large enough to 
justify the 20% loss in photons when using the r = 20" 
region. In the following we will then use the 40" radius 
spectrum, unless cxplicitely stated. 
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Table 1. Best fit results. Continua 
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4. The E>4 keV XMM-Newton spectrum 

The soft X-ray high energy resolution line spectrum of 
NGC 1068 is discussed in detail by Kinkhabwala et al. 
(2002) and Ogle et al. (2003). We limit ourselves to anal- 
yse and discuss the high energy (E >4 keV), moder- 
ate resolution EPIC-pn spectrum of the nuclear region 
of NGC 1068, with particular emphasis on the iron and 
nickel line properties. Let us, however, first briefly discuss 
the continuum properties. 

4.1. The continuum 

Following Iwasawa et al. (1997) and Matt et al. (1997a), 
we fitted the continuum with two components, describ- 
ing reflection from both neutral and ionized matter (the 
latter component including both the 'warm' and 'hot' re- 
gions, Bianchi et al. 2001). For the neutral reflection com- 
ponent we adopted, as customary, the model PEXRAV, 
with R fixed to -1 (i.e. with the primary emission shut 
off). The ionized reflector component is instead modeled 
by a simple power law, as the ionized reflector should 
act as a "mirror" for the primary emission. The power 
law index has been assumed to be the same for the two 
models. To fit the spectrum, we of course added as many 
lines as required, as discussed in details in the next para- 
graphs. The fit is fully acceptable (x 2 / d - oi -= 14 3-7/149). 
The best fit parameters for the continua, for the neu- 
tral and the ionized lines are reported in Tables [21 and 
El respectively (all errors refer to 90% confidence level 
for one interesting parameter, unless explicitely stated). 
The folded spectrum and best fit model, along with the 
model/data ratio, are shown in Fig. |3 while the unfolded 
spectrum is shown in Fig. 0] The power law index is 
2.041qq3, and the ratio between the 4-10 keV cold and 
ionized reflectors fluxes is about 1.5. (Fitting the spec- 
trum extracted from a 20" radius, a slightly steeper in- 
dex, 2.10^0 04' i s found; the ratio between the two compo- 
nents does not change significantly. All other parameters 
are also consistent with those obtained with the r=40" ex- 
traction radius, but with of course somewhat larger error 
bars). This result could be compared with that derived by 
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Fig. 3. Spectrum and best fit model, and model/data ra- 
tio. See text for details. 

Matt et al. (1997a) from the December 1996 BeppoSAX 
observation, where the flux ratio was instead 0.5, sug- 
gesting significant variations in likely both components. 
Indeed, variability on a time scale of 1 year, probably 
in the ionized reflector component, has been reported by 
Guainazzi et al. (2000). However, the BeppoSAX result 
must be taken with caution: the FLoid/^lon ratio is very 
sensitive to the power law index, which in BeppoSAX 
was only loosely determined (1.74lg'gg). We re-fitted the 
BeppoSAX spectrum imposing the power law index found 
from XMM-Newton, i.e. 2.04. The result is only slightly 
worse statistically A% 2 =2.7 with one less free parame- 
ter, corresponding to a 90% confidence level according 
to the F-test), and the ratio now becomes ~1.3. On the 
contrary, fitting the XMM-Newton spectrum imposing 
the BeppoSAX best fit power law index gives a signifi- 
cantly worst fit (x 2 /d.o.f.=174. 4/150); the ratio is now 
about 1.15. Indeed, the XMM-Newton 90% confidence 
level lower limit on T (i.e. 2.01) corresponds to a flux ratio 
of 1.49. At this confidence level, therefore, the BeppoSAX 
and XMM-Newton estimates are inconsistent with each 
other. There is therefore some, even if still not conclu- 
sive given also the uncertainties when comparing dif- 
ferent instruments, evidence for variability between the 
BeppoSAX and the XMM-JVewtori observations, taken 
3.5 years apart, which would give an upper limit to the dis- 
tance of the spatially unresolved reflecting regions of the 
order of a few parsecs. Clearly, a second XMM-A^euifon 
observation is necessary to confirm this finding. 

Another interesting result is the large iron abundance, 
Ap e =2.4 (with respect to the solar value of Anders & 
Grevesse 1989). In the fit, Ap c is basically derived from 
the depth of the iron edge in the Compton reflection com- 
ponent, which depends on the spectral index and, for 
a plane parallel geometry, also on the inclination angle. 
Solar abundances for the other elements (in particular for 
CNO elements, which contribute most to the photoelec- 
tric cross section at 6.4 keV) have been assumed. The fit 
has been performed keeping the inclination angle i fixed 
to 63°, i.e. the Xspec default value. For lower values of 
the inclination angle a less deep edge is expected, and 
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Fig. 4. Unfolded spectrum and best fit model. See text for 
details. 




Fig. 5. Ap c ~r contour plot. 



therefore an even larger value of the iron abundance is 
found (e.g. 3.3 for an inclination of about 20°). In Figs. [3] 
the Ape-T (with i free to vary) is shown; it is important 
to note that the iron abundance is larger than 1.5 at the 
99% confidence level (two interesting parameters). 

4.2. Emission lines from the cold reflector 

Of the several lines observed in the high energy spectrum 
of NGC 1068 (see Fig. |3J, four can be attributed to neu- 
tral Fe and Ni atoms, i.e. the Ka and K/3 lines from both 
elements (see Table EJ- All lines, both neutral and ion- 
ized, have been fitted with a gaussian profile with a fixed 
to 1 eV (i.e. much less than the energy resolution; no im- 
provement is found leaving a free to vary for the brigthest 
lines). All lines we attribute to the neutral reflector are 
significant at more than 99.99% confidence level, accord- 
ing to the F-test, apart from the Ni K/3 at 8.25 keV which 
is significant at the 99.90% confidence level. 

The iron Ka line energy (6.424±0.001 keV; note that 
the statistical error is much lower than the calibration un- 



certainty, which is of the order of 10 eV 1 ) is larger than the 
value for neutral iron (actually a doublet with energies of 
6.391 and 6.405 keV, House 1969, with a 1:2 branching ra- 
tio), and would instead suggest Fe xvi; however, the mea- 
sured K/3/Ka flux ratio for iron (including only the line 
core for the Ka line, see below) is 0.20^q'q3j slightly larger 
than the expected value even for neutral value, which is 
about 0.16 (see Molendi et al. 2003). This value stays al- 
most constant up to Fe ix, then starts decreasing (Kaastra 
& Mewe 1993). It is therefore unlikely that iron is much 
more ionized than that (the line energy is also basically 
constant over this ionization range, House 1969). The dis- 
crepancy with the observed line energy may result from 
a less than perfect energy calibration of the pn, or a real 
blueshift of the line, which is however hard to explain in 
the classical torus model for the neutral reflection. It is 
worth noting that the MOS2 gives exactly the same best 
fit line energy, and that Chandr a /HETG also measured 
a blueshift, albeit smaller (6.41 1±0. 007 keV, Ogle et al. 
2003). The K/3/Ka ratio for nickel is 0.57±0.35, definitely 
larger than the expected value, which is similar to that for 
iron (Kaastra & Mewe 1993). While the measured energy 
of the line, 8.25^q'q3 keV, is perfectly consistent with the 
expected value of 8.265 keV for the neutral Ni K/3 line 
(Bearden 1967), the Fe xxvi K/3 lines at energies cluster- 
ing around 8.25 2 are also possibly contributing. 

The Ni-to-Fe Ka flux ratio is 0.13±0.03, much larger 
than the expected value of 0.04-0.06 (calculated as de- 
scribed in Molendi et al. 2003 assuming a factor 2 over- 
abundance for both iron and nickel) . This implies a nickel- 
to-iron overabundance of a factor ^2, or more 

Clear evidence for an iron Ka Compton shoulder (here- 
inafter CS) is also found, confirming and refining the 
ASCA (Iwasawa et al. 1997) and Chandra/RETG (Ogle 
et al. 2003) results. The CS is fitted with a gaussian profile 
with the centroid energy and a fixed to 6.3 keV and 40 
eV, respectively (Matt 2002). If the CS is not included a 
significantly worse fit (x 2 /d.o.f.=181. 2/149) is found, even 
allowing for the width of the main line to vary. The ratio 
between the CS and the line core fluxes is about 0.2, in 
agreement with the theoretical expectation for Compton- 
thick reflecting material (Matt 2002) , so suggesting an ori- 
gin in the visible part of the Njj >10 25 cm~ 2 absorbing 
material (Matt et al. 1997a). 

The iron Ka EW (with respect to the cold reflection 
continuum) is 1.2 keV (1.45 keV including the CS), rea- 
sonably in agreement with theoretical expectations (e.g. 
Ghisellini et al. 1994, Matt et al. 2003) even with the 
observed iron overabundance, given the less than linear 
increase of the line EW with the iron abundance (Matt et 
al. 1997b). 



see http: //xmm. vilspa. esa. es/docs/do cuments/CAL-TN-00 18-2-1 .pdf 

2 Unless explicitely stated, line energies are 
taken from the NIST Atomic Spectra Database, 
http: //physics .nist . gov/cgi-bin/AtData/main_asd. 
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Table 2. Best fit results. Lines from cold matter. 
Equivalent widths are calculated against the cold reflector 
only. 



4.3. Emission lines from the ionized reflector 

Four emission lines from highly ionized matter (the 'hot' 
reflector cited above) are also detected in the pn spectrum 
(Table[2J). All of them are significant at more than 99.99% 
confidence level, according to the F-test, apart from the Ni 
He-like Ka at 7.83 keV which is significant at the 99.79% 
confidence level. Besides the He- and H-like Fe Ka lines 3 , 
already discovered by ASCA (Ueno et al. 1994; Iwasawa et 
al. 1997; Bianchi et al. 2001), two lines at about 6.61 and 
7.83 keV are also detected with high significance. The en- 
ergy of the first line points towards fluorescence from Be- 
like iron, which however has a quite low fluorescent yield 
of 0.11 (Krolik & Kallman 1987). Moreover, the matter 
may likely be optically thick to the line, which results in 
an effective Auger destruction. It is therefore possible that 
this line is actually a blend of Be- and Li-like lines (the 
latter having an energy of about 6.65, a high fluorescent 
yield, 0.75, and no possibility of Auger destruction). (The 
possibility that this line is the CS of the He-like iron line 

3 Also for these lines the measured energies are not in agree- 
ment with the expected values, if only statistical errors are 
taken into account. However, in these cases identification is 
certain, as no other strong lines are present at those energies. 
It is worth noting that the energy of the He-like line is greater 
than expected, and that of the H-like line lower, ruling out 
inflow or outflow if the two lines, as likely, are produced in the 
same region. 



is ruled out by the low optical depth of the hot reflecting 
matter, as discussed below). The line at 7.83 keV is likely 
a blend of the Ni He-like Ka (7.81 keV) and Fe He-like 
K/3 lines (7.87 keV). The rather wide range of ionization 
argues against collisional ionization. On the contrary, we 
have verified with CLOUDY 4 that photoionization may pro- 
duce such a range in a single zone with e.g. a column den- 
sity a few times 10 21 cm -2 (see below) and a ionization 
parameter of about 150. 

The measured values of the EWs can be compared with 
the results of the best previous observation, i.e. the 100 
ks ASCA observation (Bianchi et al. 2001). The results 
reported there are somewhat different, but were obtained 
with a simpler model. We refitted the ASCA spectrum 
with all the lines found by XMM-TVeuifon (and the power 
law index fixed to 2.04) and found that, within the errors, 
the fluxes measured by ASCA are consistent with those 
measured by XMM-Newton. It should also be noted that 
the presence of the 6.61 keV line vindicates the claim of 
Iwasawa et al. (1997), based on the first ASCA observa- 
tion, of a redshift of the ionized lines. 

The EW of the He- and H-like iron lines can be com- 
pared with the results of Matt et al. (1996) to derive at 
least an order-of-magnitude value for the column den- 
sity of the reflecting material. In their figures 5 and 6 the 
EW of these lines, calculated including contributions from 
both recombination and resonant scattering, are shown as 
a function of the column density of the reflecting mate- 
rial. The calculations assumed iron abundance and frac- 
tions equal to 1 of the considered ions. The measured iron 
overabundance roughly compensates for the less than one 
fraction of each ion (for instance, our CLOUDY calculation 
gives a value of 0.4 for the fraction of He-like iron atoms, 
and somewhat smaller value for the H-like atoms). Even 
taking into account the fact that the ionized continuum 
includes also reflection from the warm reflector (Bianchi 
et al. 2001 estimated that the hot reflector contributes to 
about 3/4 of the total ionized reflection), from the above- 
mentioned figures of Matt et al. (1996) a column density 
of the hot ionized reflector of a fewxlO 21 cm -2 can be 
deduced. 



5. Summary 

We have analysed the high energy (>4 keV) XMM- 
Newton EPIC-pn spectrum of NGC 1068. The main re- 
sults can be summarized as follows: 

i) Possible (but not conclusive) evidence for variations 
of both neutral and ionized reflectors with respect to the 
December 1996 BeppoSAX observation is found, suggest- 
ing an upper limit to the size of both reflectors of a few 
parsecs. Another XMM-Newton observation, to compare 
results obtained with similar instruments, is however re- 
quired to confirm this finding. It is worth noting that a 
variability on a time scale of about a year of the ionized 

4 http://www.nublado.org/ 
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Table 3. Best fit results. Lines from ionized matter. 
Equivalent widths are calculated against the ionized re- 
flector only. 



reflector was found by Guainazzi et al. (2000) comparing 
two BeppoSAX observations. 

ii) Iron is overabundant, with respect to lower Z ele- 
ments and when compared to the solar value, by a factor 
about 2. Nickel is, in its turn, overabundant by a factor 
~2 with respect to iron. A qualitatively similar result was 
found in the Circinus Galaxy (Molendi et al. 2003), where 
however both overabundances were smaller. 

iii) The iron K/3/Ka flux ratio points to a low ion- 
ization state of iron, inconsistent with the Ka line en- 
ergy, suggesting calibration problems. The nickel K/3/Ka 
is larger than expected, but the K/3 line may actually be 
a blend with ionized iron lines. 

iv) The Fe Ka Compton Shoulder is also detected, with 
a relative flux of 0.2, in agreement with the value expected 
for reflection from Compton thick matter (Matt 2002) 

v) Be-, Li-, He- and H-like iron emission lines, as well 
as He-like Ni line, are also found. Their EW suggests a 
column density of the hot ionized reflector of a fewxlO 21 
cm~ 2 (Matt et al. 1996). 
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